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Abstract. The Monte Carlo HYDJET++ model, that contains both 
hydrodynamic state and jets, is applied to study the influence of the interplay 
between soft and hard processes on the formation of the elliptic flow in heavy- 
ion collisions at RHIC and LHC energies. Jets are found to cease the hydro-like 
increase of the elliptic flow with rising pt- Since jets are more influential at LHC 
than at RHIC, the V2 at LHC should be weaker than that at RHIC. Violation 
of the number-of-constituent-quark (NCQ) scaling is predicted. The decays of 
resonances are found to enhance the low-p^ part of the elliptic flow of pions and 
light baryons, and work toward the fulfillment of the NCQ scaling. 
PACS numbers: 25.75.-q, 25.75.Ld, 24.10.Nz, 25.75.Dw 



1. Introduction 

The elliptic flow of hadrons, produced in the course of non-central heavy-ion collisions, 
is considered as one of the main probes of the reaction early stage. It originates from 
the rescattering process in the initially anisotropic overlapping region. Recall that 
elliptic flow is the second coefficient of Fourier decomposition of particle invariant 
cross section in the azimuthal plane, namely v% — (cos (20)) = ((p^ — pfy/pj)- Here 
<p, pt, Pxj and p y are the azimuthal angle, particle transverse momentum, and its 
in-plane and out-of-plane components, respectively. The spatial anisotropy of the 
almond-shaped area is then transformed to the anisotropy in the momentum space. 
When the isotropy in the coordinate space is restored, the elliptic flow stops to develop 
anymore. Thus V2 is linked to the equation of state (EOS) of hot and dense partonic 
matter formed in the very beginning of the collision. 

Many features of the flow at energies up to top RHIC energy yfs = 200 AGeV 
are successively reproduced by several models, in particular, hydrodynamic models 
(for review, see (T] and references therein). However, predictions of these models 
for the elliptic flow at top LHC energy y/s = 5.5ATeV often contradict each 
other. For instance, simple scaling model [2] favors increase of vi at midrapidity 
to vi(y = 0) sa 7 — 8%, several models predict saturation of the flow [2], whereas ideal 
hydro indicates [3j|4] that py-integrated flow at LHC will be larger than that at RHIC 
although V2 HC (pt) < vf mc (pr) in the domain of low and intermediate transverse 
momenta. 
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Jets contribute to the elliptic flow via the so-called jet quenching mechanism that 
accounts for the quite weak flow of about 5%. The simultaneous description of the 
flow caused by the soft processes and the jets, and study of the interplay between 
these two mechanisms has been done recently in [5] by means of the HYDJET++ 
model [6 . The model contains the soft physics part, represented by the parametrized 
hydrodynamics with given freeze-out conditions [7], and the jets. The production of 
hard multiparton states and their propagation in dense partonic matter takes into 
account radiative and collisional losses [8] associated with the parton rescattering 
in the expanding QGP. The contribution of both soft and hard parts to the total 
multiplicity of a heavy-ion collision at ultrarelativistic energies depends merely on 
the collision energy and centrality. The few free parameters of the model have been 
tuned to provide adequate description of hadron spectra, elliptic flow, femtoscopic 
momentum correlations, and high px tails of hadron spectra in gold-gold collisions at 
RHIC energies. Further description of the HYDJET++ model can be found in [BHH]. 
The main aim of our present paper is the study of the interplay between soft and 
hard processes responsible for the elliptic flow production in heavy-ion collisions at 
energies of RHIC and LHC, respectively. The influence of the resonance decays on the 
final elliptic flow of particles is very important at both energies in question. Here, we 
are going to benefit from the rich table of baryon and meson resonances (about 360 
states) implemented in the HYDJET++ model. 

2. The influence of jets and resonances on the V2 formation 

The model description of the particle elliptic flow as a function of the transverse 
momentum is plotted in Fig. [1] onto the RHIC experimental data taken at y/s = 
200AGeV. Besides the standard rise of the flow excitation functions and the mass 
ordering of the particle flow at < pr < 2 GeV/c, the model exhibits also the 
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Figure 1. The V2 (pt) dependence in the HYDJET++ model for different hadron 
species (lines) and comparison with RHIC data (symbols). 

saturation and the drop of the particle ^(pt) at pt > 2.5GeV/c accompanied by 
the violation of the mass ordering. The last two features cannot be attributed to 
ideal hydrodynamics that indicates the instant rise of the elliptic flow up to unity 
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without the crossing of meson and baryon branches as pt rises to infinity. Here we 
see the result of the interplay between the soft hydro processes and hard jets. In 
contrast to hard processes, the multiplicity of hadrons produced in soft processes 
drops exponentially with increasing transverse momentum. Therefore, after a certain 
Pt the particle spectrum is dominated by the jet particles. The latter carry quite weak 
elliptic flow, thus the combined flow of high pr hadrons drops. Also, the slope of the 
pr spectra of heavy hadrons is not as steep as that of the light particles. Because of 
this, the hydro component of the transverse momentum distribution of heavy particles 
dominates until larger values of pt resulting to the change of the mass ordering of the 
flow: After a certain pt the heaviest particles possess the largest flow. 

What are the consequences of such peculiarity in the behavior of elliptic flow 
for the heavy-ion collisions at LHC energies? Will the hadronic elliptic flow extend 
its strength or not? To answer these questions we generated 1,000,000 Au + Au 
collisions at y/s = 200 AGeV and 500,000 Pb + Pb collisions at y/s = 5.5 ATeV, both 
reactions at fixed centrality c = 42%. The centrality was chosen to observe already 
a strong elliptic flow and, simultaneously, to keep the flow fluctuations at a rather 
modest level. Transverse momentum dependence of elliptic flow of charged pions, 
(anti)protons, charged kaons, and A + S is displayed in Fig. [2] for Au + Au and Pb + 
Pb collisions at RHIC (left panels) and LHC (right panels) energies, respectively. One 
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Figure 2. Left: The V2(pt) dependences of all hadrons (full symbols) and of 
direct hadrons (open symbols) in the HYDJET++ simulations of Au + Au at 
v/J = 200 AGeV with centrality c = 42%: (a) n ch , (b) p + p, (c) K ch , and (d) 
A + S. Right: Similar distributions, but for Pb + Pb collisions at y/s = 5.5 ATeV; 
(e),(f),(g),(h) correspond to (a),(b),(c) and (d), respectively. 



can sec that the flow of the most abundant hadron species at RHIC is stronger than 
its LHC counterpart. Although the soft hydro parts of the elliptic flow at RHIC and 
at LHC are similar [5], the hard parts are not. More particles are produced at LHC 
via the jet fragmentation. Therefore, hadrons from jets start to dominate over the 
hadrons, originated from the soft processes, at lower pt as compared with RHIC. The 
resulting elliptic flow at LHC is thus severely reduced. Another fact distinctly seen in 
Fig. [2] is the influence of the resonance decays on the elliptic flow of directly produced 
particles. Note that according to the model calculations, resonances contribute to the 
production of 80% of pions, 70% of protons, 60% of A + E, and more than 50% of 
kaons at LHC. For all hadrons, except kaons at LHC, final state interactions (FSI) 
increase the particle flow in the interval 1.5GeV/c <pt< 6GeV/c. However, decays 
of resonances diminish the pion flow at pt > 1.5 GeV/c for both energies in question. 
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The contributions of A decays to the flows of pions and protons, and p and ui 
decays to pion elliptic flow are presented in Fig. [3] Elliptic flows of all resonances are 
weaker than that of pions at low transverse momenta and stronger at px > 1.5 GeV/c. 
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Figure 3. Left: The i>2(pt) dependence for (a) ir ch and (b) p + p produced both 
directly (circles) and in A decays (squares) in the HYDJETH — h simulations of 
Au + Au at v 7 ! = 200AGeV with centrality c = 42%. (c) and (d) The same as 
(a) but for charged pions produced in decays of p and ui, respectively. The flow 
of resonances is shown by triangles. Right: Similar distributions, but for Pb + 
Pb collisions at y/s = 5.5ATeV; (e),(f),(g),(h) correspond to (a),(b),(c) and (d), 
respectively. 



Due to the decay kinematics, pions from the decay uj — > 3n are getting softer 
Pt distributions compared with that of to. In contrast, the transverse momentum 
distributions of pions from the p — > 2ir decay are close to that of the p meson. As 
a result, elliptic flow of pions from the uj decay is transferred to the softer pr region 
compared with v% (pt), whereas elliptic flow of pions from the p decay is similar to 
v 2(pt)- Decays of A's produce protons with transverse momenta close to that of the 
mother particles and pions with rather soft pr spectra. Again, the elliptic flow of 
pions is shifted to softer px region. Therefore, the resulting elliptic flow of pions is 
a bit lower than the vi of directly produced pions at pt > 1.5GeV/c and stronger 
at intermediate transverse momenta. Modification of the elliptic flow of pions and 
light baryons can lead to violation of the hydro-induced mass hierarchy in the V2{pt) 
sector. Particularly, the so-called number-of-constituent-quark (NCQ) scaling [9l[T0] 
can be violated. Below we check this hypothesis. 



3. Violation of NCQ scaling 



The NCQ scaling claims the similarity of the V2 / n q (K Et / n q ) distributions, where 
both the elliptic flow V2 and the transverse kinetic energy of a hadron KEt are divided 
by the number of constituent quarks n q . This scaling holds up until KET/n q w 1 GeV 
[TTj . The HYDJET++ calculations of the V2/n q (KET /n q ) in heavy-ion collisions at 
RHIC and LHC are shown in Fig. 0] both for directly produced particles and for 
all hadrons. One can see that for direct hadrons the NCQ scaling is approximately 
satisfied at KET/n q < 1 GeV at RHIC only. At LHC the NCQ scaling is not achieved 
for the direct particles due to the strong influence of the jets. To see possible deviations 
more distinctly, the particle distributions are also normalized in Fig. [4] to the elliptic 
flow of lambdas, v^/n 1 * : 
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Figure 4. Left: Upper row: The KEx/nq dependence of elliptic flow for (a) 
direct hadrons and (b) all hadrons produced in the HYDJETH — h simulations 
of Au + Au at = 200AGeV with centrality c = 42%. Bottom row: The 
KEx/riq dependence of the ratios (v2/n q ) j /3) for (c) direct hadrons and 
(d) all hadrons. Right: Similar distributions but for Pb + Pb collisions at 
•ys = 5.5ATeV; (e),(f),(g),(h) correspond to (a),(b),(c) and (d), respectively. 



The decays of resonances essentially improve the situation. For instance, the 
elliptic flow of mesons has no resonance feed-down, whereas the flows of protons and 
lambdas are significantly enhanced by heavy resonances at KEt > 0.5 GeV. Because 
of the resonance decays, the NCQ scaling holds at 10% up to KEt > 1.0 GeV for 
RHIC energies. At LHC the effect of jets is too strong, and only the approximate 
number-of-constituent-quark scaling is realized. 

4. Conclusions 

In conclusion, the interplay between hydrodynamic processes and jets, and its influence 
on the formation of the particle elliptic flow is studied within the HYDJET++ model 
for heavy-ion collisions at RHIC and LHC energies. Jet particles that carry quite 
weak elliptic flow start to dominate over the hydrolike particles at a certain transverse 
momentum, thus causing the reduction of the combined i>2. Therefore, we predict 
that elliptic flow at LHC will be smaller than that at RHIC, v^ HC (p T ) < vf mc (p T ), 
at p T > 3 GeV/c. 

Besides, jets also account for reversing of the mass ordering of the elliptic flow 
at intermediate and high transverse momenta: Here the heaviest hadrons possess the 
largest flow. Finally, the decays of resonances are shown to push the hadron excitation 
functions V2/n q (KET /n q ) toward the fulfillment of the NCQ scaling. However, elliptic 
flow of jet particles at the LHC energy should interfere with its hydrolike counterpart 
already at intermediate transverse momenta. Thus, the realization of the approximate 
number-of-constituent-quark scaling becomes worse compared with the RHIC case. 
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